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The Kondo effect in a + ipj^-wave superconductor is studied by applying the Wilson's nu- 
merical renormalization group method. In this type of superconductor with a full energy gap 
like a s-wave one, the ground state is always a spin doublet, while a local spin is shrunk by 
the Kondo effect. The calculated magnetic susceptibility indicates that the spin of the ground 
state is generated by the orbital effect of the px + ipy-wave Cooper pairs. The effect of spin 
polarization of the triplet superconductor is also discussed. 

KEYWORDS: Kondo effect, unconventional superconductivity, triplet Cooper pair, numerical renormalization 
group 



§1. Introduction 

The Kondo effect is strongly suppressed by existence 
of an energy gap at the Fermi energy, since infinitesimal 
excitations are important. In fact, two energy scales, the 
Kondo temperature Tk and the energy gap A, compete 
with each other. In a standard BCS (s-wave) supercon- 
ductor, the competition is clearly seen in Tk/A d( 
dence of a bound state energy below the energy gap.El 
The ground state is a spin doublet for Tk/A <C 1 and 
it is a spin singlet for Tk/A ^ 1. The interchange of 
ground states occurs at around Tk/A = 0.3. 

According to the recent nuclear magnetic resonance j-, 
measurement performed in Li-substituted YBa2Cu307_5,El' 
the local spin susceptibility indicates the reduction of 
the Kondo effect due to opening of the superconducting 
energy gap below the superconducting transition tem- 
perature. This is supported by the recent theoretical 
study an the Kondo effect in a d^2_y2-waNe supercon- 
ductorJa^ It seems that the interchange of spin-doublet 
and spin-singlet ground states also occurs in such singlet 
superconductors as the s-wave and d-wave. 

In this paper, we study the Kondo effect in triplet 
superconductors, using- the mmierical renormalization 
group (NRG) methodB'LP The simplest form for the 
NRG calculation is & px -\- ip^-wave pairing type with 
a full energy gap like the s-wave. One of the Px + ipy- 
wave pairing states is given by an odd vectorial function 
d = z{kx + ifcy) or is expressed more^fixplicitly with a 
gap function A^a'ik) = A||(fca; +ifcy).El' A possibility of 
this state has been discussed for a triplet superconduc- 
tor Sr^Ru04 in which the time-reversal symmetry is bro- 
ken.Qll3^ We also study the case of a spin-polarized super- 
conducting state represented by both A||(fcj; -I- iky) and 
/S.ii{kx+\ky). This state is a possible candidate proposed 
for a ferromagnetic superconductor recently discovered 
in UGe2 ,liil'li3' where such a triplet superconductivity is 
expected to be realized under high pressure environment. 
This is also expected in other ferromagnetic metals such 
as ZrZnzEJ) and URhGclli 



Since the Cooper pair of the px + ipy-wave state has a 
net angular momentum Z = 1, it is formed by conduction 
electrons having I — and I — 1 angular momenta. A 
local spin is coupled directly to the I = conduction 
electrons. When the superconducting state is set on, the 
/ = 1 conduction electrons also participate in the Kondo 
effect. In this case, tha^ground state is a spin doublet 
for all values of Tk/AJ13^ which is different from the s- 
wave case. In order to elucidate magnetic properties of 
this new type of spin-doublet ground state, we calculate 
magnetic susceptibilities of the impurity. In the Kondo 
problems, two kinds of magnetic susceptibilities are often 
discussed: one is a local spin susceptibility xiocai and the 
other is called an impurity susceptibility Ximp- The size 
of the local spin is reduced by the Kondo effect as Tk/A 
increases, which is reflected in the zero temperature limit 
of Txiocai- In such gapped systems as superconductors, 
it is qualitatively different from Tximp- The difference is 
also found in semi-conducting, systems where a pseudo- 
gap exists at the Fermi energy.llaO'EJ) It is more evident 
in the fully gapped system we discuss here. 

We also discuss the spin effect of Cooper pairs in the 
triplet superconductors by introducing spin polarization 
(|Ai--|-| — |Aj^||). When either A|| or A|| vanishes, the 
infinitesimal excitations are the most relevant and the 
Kondo effect becomes inactive. In this case, the spin 
direction of low-energy excitations is locked by the zero 
order parameter, and the local spin behaves as a classical 
spin. 

This paper is organized as follows. In §2, we present 
the Kondo model for the px + ipy superconductivity and 
the corresponding NRG Hamiltonian. In §3, we show 
the NRG results for spin-unpolarized and spin-polarized 
superconducting states separately. Concluding remarks 
are given in §4. 

§2. Model Hamiltonian 

Let us begin with the following Hamiltonian for 
the spin-unpolarized triplet Cooper pair represented by 
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(2.4) 
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Here akia is an annihilation operator of the conduction 
electrons, where the subscripts k, I, and a represent 
wave number renormalized by fcp, angular momentum 
and spin, respectively. J (< 0) is an antiferromagnetic 
exchange coupling constant. The dispersion relation of 
the conduction electrons has been linearized for each an- 
gular momentum / in eq. (2.2). The total angular mo- 
mentum of the -I- ipy-wave Cooper pair is equal to one, 
which is given by the coupling of two angular momenta, 
I and —l+l {I — 0, ±1, ±2, • • •), as shown in eq. (2.3). 
The local spin is coupled directly to the I — electrons. 
Therefore the relevant Px + ipj,-wave Cooper pair in the 
Kondo effect consists of the I — and I = 1 conduc- 
tion electrons. In the normal metallic state (A = 0), the 
I = conduction electrons screen the local spin, while 
the I = 1 conduction electrons are completely decou- 
pled from the local spin. In the superconducting state 
(A 0), the I = 1 conduction electrons participate in 
the Kondo effect. 

As we have derived in our previous paper the NRG 
Hamiltonian for the Px+ipy-yf&ve is given in the following 
recursion relation: 



Hn+1 = K^'^Hn + ^ \eN{.C^N+l.ra'^Nra + H.C.) 
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J and A are the effective exchange coupling J and the su- 
perconducting energy gap A normalized by (1-1- A~"'^)/2. 
A is a logarithmic discretization parameter and it is taken 
to be A = 3 throughout this paper. CNra is an operator 
of the NRG fermion quasiparticle in the A^-th shell. The 
subscript r = ± represents the two channels which are 
constructed by the / = and I = 1 orbitals. The Kondo 
temperature is defined by Tk = | J|^^^ exp(— 1/| J|). 

In the same manner, we can derive the NRG Hamilto- 
nian for describing the Kondo effect in the spin-polarized 
Px + ipy-wave state represented by A|| and A||. Since 
the BCS interaction corresponding to eq. (2.3) is given 
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Fig. 1. Tk/A dependence of the bound state energy (_Eb/A)* 
and the effective g-factor for the spin-unpolarized superconduc- 
tivity, go is a bare g-factor of the local spin. The lowest-lying 
1200 states are kept at each renormalization step in the NRG 
calculation. 
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we only put spin dependence into the energy gap as A ^ 
A^a- in eqs. (2.6) and (2.7). 

Next we describe how to calculate magnetic suscep- 
tibilities. There are two kinds of them. One is a local 
impurity susceptibility defined by 



Xlocal(r) = gflB- 



Tr [g,exp(~/3gjv)] 
Tr [cxp(-/3i?Ar)] 



(2.9) 
(2.10) 



Here Sz is the local spin operator, and /3 ~ 1 is 
taken to reduce cutoff dependence in the NRG calcu- 
lation. The temperature T is defined by /3T = [{I + 
A-1)/2]A-(^-i'/2Me|) a small magnetic field h is ap- 
plied only at the local spin site. In practice, the impurity 
part of the Hamiltonian Hq in eq. (2.7) is replaced by 
Hq — g^B'S'z^A^^/^ in the NRG calculation. The other 
susceptibility is defined by 



Ximp(T) = x{T) - Xo{T), 



(2.11) 



where x ixo) is a magnetic susceptibility of the conduc- 
tion electrons including (in the absence of) the local spin. 

§3. Results 

In the remaining parts of this paper, we show the 
NRG results for the spin-unpolarized (A|| = A||, A|| — 
All = 0) and spin-polarized (A|| ^ A^, A|| = A^^ = 
0) superconducting states. 



3.1 Spin-unpolarized suimfconductivity 

In our previous study,E£r we found that the ground 
state is always a spin doublet {S — 1/2) in all the 
Tk/A region. The first excited state (bound state) is 
a particle- hole doublet with no spin (5 = 0). The ratio 
of the bound state energy and the superconducting en- 
ergy gap converges after the renormalization process. We 
express the convergent value with {Eb/ A)* . It is scaled 
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Fig. 2. Temperature dependence of Txiocal in th^ strong coupling 
region. The unit is (go/^B)^! where go is a bare g- factor of the 
local spin. The data for Tk/A > 8 are plotted here: the circles 
(A = 1.5A = 5 X 10~^), squares (1 X 10~*), and diamonds 
(2 X 10-*) are for the fixed Tk = 1.76 X IQ-^; the up triangle 
(A = 2 X 10"^), left triangle (5 X 10"^), and down triangle 
(1 X 10""') symbols are for the fixed Tk = 5.34 X 10~*. In the 
middle region, Xlocal shows a constant value (oc T^^). The zero 
temperature limit gives the effective g-factor (oc A/Tk) in Fig. 1. 



by Tk/A, as shown in Fig 1. As Tk/A increases, the 
bound state {S = 0) energy approaches that of the spin- 
doublet ground state {S = 1/2). It implies that the local 
spin is quenched at Tk/A — > oo. The screening of the 
local spin is clearly found in the zero temperature limit 
of Txiocab exhibiting the Curie law. Therefore this value 
corresponds to an effective g-factor of the local spin. In 
Fig. 1, the bound state energy and the g-factor exhibit 
similar Tk/A dependence. In the strong coupling region 
(Tk/A » 1), both of them behave as A/Tk- In addi- 
tion, the temperature dependence of Txiocai for T < Tk 
is scaled nicely as shown in Fig. 2. This can be expressed 
by the following equation: 



Txiocal 

(5oMb)^ 



A 



/ 



TkT 
A2 



(3.1) 



where go is a bare (/-factor of the local spin and the func- 
tion f{x) behaves as 



f{x)^x (1< X < (Tk/A)2), 
f{x) ~ 10^1 (0 < X < 1). 



(3.2) 
(3.3) 



The constant xiocai ^ in the middle region corre- 
sponds to the low temperature behavior in the metallic 
case (A = 0). On the other hand, at high temperatures 
(T > Tk), Txiocai is scaled by the single energy Tk as 
expected by the usual Kondo effect, and the energy gap 
A can be neglected. As Tk/A decreases, the middle re- 
gion of Txiocai is narrower. For the weak coupling limit 
(Tk/A <^ 1), Txiocai is not reduced so much at high tem- 
peratures (T > A), and it becomes constant near to 1/4 
for T < A, corresponding to the value for a free S" = 1/2 
local spin. Thus the competition between the Kondo ef- 
fect and the energy gap is observed clearly by the zero 
temperature limit of Txiocai (the effective g-factor) and 
its temperature dependence for Tk/A ^ 1. 

Next we discuss the two kinds of magnetic susceptibil- 
ity Xlocal and Ximp. They show the different temperature 




Fig. 3. Temperature dependence of Tximp and Txiocal ior the 
spin-unpolarized superconductivity. Tk and A{= 1.5A) are cho- 
sen as Tk = 5.34 x 10"'' and A = 10~^, respectively. The 
lowest-lying ^ 1200 states are kept at each renormalization step 
in the NRG calculation. 



dependence as seen in Fig. 3. At higher temperatures, 
we can see that both exhibit similar temperature depen- 
dence. For Tximp, it upturns at sufficiently low tem- 
peratures and deviates from Txiocai- It saturates at low 
temperatures and takes 1/4 which is the same value with 
Txiocal for a free S' = 1/2 local spin. On the other hand, 
Txiocal in Fig. 3 is reduced by the Kondo effect. What 
does this difference mean? Let us discuss this point in the 
strong coupling limit (Tk/A ^ oo) where Txiocai — > 
at low temperatures. Although the ground state keeps 
a spin doublet, the effective g-factor of the local spin 
is strongly reduced by the Kondo effect (see Fig. 1). 
This implies that the I — conduction electrons quench 
the local spin by forming a Kondo singlet. In order to 
gain the superconducting condensation energy, the I = 1 
conduction electrons form Cooper pairs with the I — 
electrons. As a result, the I = 1 electrons are connected 
with the Kondo singlet. This concludes that the spin 
of the ground state is generated by the I = 1 conduc- 
tion electrons in the strong coupling limit. For xiocai, 
the I = 1 conduction electrons do not respond to the 
magnetic field directly, since it is applied only at the im- 
purity. On the other hand, their magnetic response is 
reflected in Ximp- Thus Ximp represents the magnetic re- 
sponse to an external field applied not only at the local 
impurity but also to the conduction electrons. The qual- 
itative difference between xiocai and Ximp is obvious in 
the gapped systems. 

In Fig. 3, we note that the upturn in Tximp takes place 
at around T = A. For T > A, the energy gap does not 
affect Ximp, so that Tximp has almost the same temper- 
ature dependence with Txiocai- Therefore we conclude 
that the different behavior between xiocai and Ximp for 
T < A is due to the finite energy gap in the density 
of conduction electron states. This should be compared 
with £L4)S£Li|Uio-gap case where such a difference is also 
found. ESEZHia^ We note that infinitesimal excitations ex- 
ist in this case, while we treat the full gap systems. If 
the density of states at the Fermi energy (e = 0) is given 
by p cx |er, both Txiocai(T ^ 0) and Tximp(T ^ 0) 
depend on r strongly. The former is extremely small, 
while the latter takes finite values between and 1/4. In 
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Fig. 4. (a) Tk dependence of the bound state energy (Eb/A)* = 
(Eb/Aj^j^)* and (b) Tk dependence of the effective p-factor. go 
is a bare value of the g-factor of the local spin. In both figures, 
All = 0.01 is taken. The circle, star, triangle and plus symbols 
represent the results for fixed A(= 1.5A) = 0.01, 0.008, 0.001 
and 0, respectively. The ground state is always a spin doublet 
and the energies of the bound state are measured from those 
of the ground state. The first excited state is a particle-hole 
doublet (5 = 0). The lowest-lying ~ 500 states are kept at each 
renormalization step in the NRG calculation here. 



general, the difference of the two susceptibiUties is due to 

the existence of an energy gap. Our study on the full-gap 
superconducting state shows this fact more evidently. 

3.2 Spin-polarized superconductivity 

The spin-polarized superconducting state is represented 
by the spin dependent order parameters A|| and An. 
For A|| = A||, we obtain the same results as in the 
unpolarized spin case. The polarization (|A||| — |A|J) 
suppress the strength of the Kondo effect, so that it in- 
creases the bound state {S = 0) energy. In Fig. 4(a), the 
effect of the difference between A|| and A^j^ is found 
as the increase of the bound state energy for each val- 
ues of Tk when Tk is smaller than A||. As expected, 
the Tk dependence is close to that for the unpolarized 
case (A|| = A^) if Tk > max(|A|||, lA^J). We also 
calculate xiocai for each value of A^^, keeping A|| fixed 
at 0.01. As shown in Fig. 4 (b), the effective fl'-factor 
decreases with the increase of Tk- It becomes smaller 
for each Tk as An is decreased from A||. As long as 
we see Fig. 4(b), the data for A^^ = 0.001 (10% of A^) 
is not much different from those for A^^ = where in- 
finitesimal excitations arc dominant. Therefore the spin 
polarization of the energy gap just modifies the effective 
5f-factor even if one of the spin dependent order param- 
eters vanishes. The local spin is completely quenched 
only when both A|| and A^^ vanish. 

The All — state mentioned above has infinite low- 
energy excitations below the energy gap associated with 
only down spin of the conduction electrons. This com- 
pletely polarized state is described by the following NRG 
Hamiltonian: 

Ht, = H^{J = 0) - J*A(^-^)/Vo\M^- 



+\[l + {-ir]S., (3.4) 

where J* is an effective coupling constant and h* is an 
effective magnetic field, which are determined by the low- 
lying excitation energies. The latter is an effect of the 
Cooper pairs with only up spin. The NRG Hamiltonian 
Hn{J = 0) for the free electrons can be described with 
quasiparticle and quasihole excitations as 

(Af+l)/2 

Hn{J = 0)= J2 ViigUi + hlhi) (AT: odd) ,(3.5) 

i=l 

Hn{J = 0) = r]'oglgo 

N/2 

+ T.'^ii9l9i + hlh^) (A^ieven), (3.6) 

i=l 

where gi and hi are annihilation operators for the par- 
ticles and holes, respectively. The energies for A = 3 
are given by 7?i = 0.800, 772 = 2.997, ■■■,Vi = 3'"\ and 
V'o = 0> V'l = 1-696, T]'^ = 5.196, = 3*-i/2_ x^e 

fermion operator at the impurity site is given by 

(Ar+l)/2 

/oi = 3-(^-i)/^ E <9i + hl) (AT: odd), (3.7) 

i=l 

f,l = 3-(^-^)/^K5o 

N/2 

+ Y.a',{gi + hl)] (AT: even), (3.8) 

i=l 

where the coefficients for A = 3 are oi = 0.628, a2 = 
0.895, ■■■,a^= Q!3('-i)/2 (q, ^ 0.5155), and a'g = 0.620, 
a[ = 0.706, •••,«■= a'3(*-i)/2 ^ 0.6784). Then 
the effective NRG Hamiltonian for describing the case of 
All = is finally expressed with 

(JV-|-l)/2 
i=l 

-r c^UaUi - h%)s, + o{[r]% (3.9) 

i 

for odd numbers of the renormalization step N and 

N/2 

H*N = Vo9o9o + Y ^i(9i9i + hjhi) 

i=l 

-ria'ofglgoS:, - J* - hlhi)S, 

i 

+h*S, + 0{[J*f), (3.10) 

for even numbers. In order to reproduce the low-energy 
excitations, J* and h* in eq. (3.10) must satisfy the fol- 
lowing relation: 

r,'n--^Hf + \--\, (3.11) 

which gives a twofold degenerate ground state (|— ) and 
g\\+)), where |±) represents the Sz = ±1/2 local spin 
state. Due to the effective field h* = (J*/2)(ao)^, all 
the NRG low-energy excitations are doubly degenerate 
as obtained for odd numbers of N. 
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If the triplet superconducting state is realized, the 

spin-polarized excitations are dominant below the higher 
transition temperature associated with a larger one of 
|A-|-|-| and |A|||. In the temperature region between the 
two transition temperatures, the local spin behaves like a 
classical spin and it is coupled to the spin-polarized con- 
duction electrons. When the second transition occurs at 
the lower temperatures, the energy gap opens completely 
and the Kondo effect generates a bound state. 

§4. Conclusion 

This study shows a variety of Kondo effect influenced 
by the triplet Cooper pairs. The Kondo effect is affected 
not only by the energy gap but also by the orbital de- 
grees of freedom of the Cooper pairs, which is completely 
different from the s-wave case. Due to the spin doublet 
ground state, xiocai displays the Curie law at suiEciently 
low temperatures, and the effective (^-factor of the lo- 
cal spin is strongly reduced in the large Tk/A region. 
The Curie term Tximp = 1/4 indicates that the spin 
of the ground state is generated by the / = 1 conduc- 
tion electrons via the orbital effect of the px + ipy-wnvc 
Cooper pairs. In addition, the spin-polarization effect is 
unique to such triplet superconducting states. If either 
or All vanishes, the local spin is not quenched by 
the Kondo effect and behaves like a classical spin. 

All of the results are useful for explaining the coexis- 
tence of magnetism and such unconventional supercon- 
ductivity. This argument is valid for heavy fermion su- 
perconductors realizing Tk/A ^ 1, while in the same 
condition the local moment vanishes at low temperatures 
for the s-wave superconductivity. 
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